Radiation from the LTB black hole 
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Does a dynamical black hole embedded in a cosmological FRW background emit the Hawking 
radiation where a globally defined event horizon does not exist? What are the differences to the 
Schwarzschild black hole? What about the first law of black hole mechanics? We face these questions 
using the LTB cosmological black hole model recently published. Using the Hamilton- Jacobi- and 
radial null geodesic- methods suitable for dynamical cases, we show that it is the apparent horizon 
which contributes to the Hawking radiation and not the event horizon. The Hawking temperature 
is calculated using the two different methods giving the same result. The first law of LTB black 
hole dynamics and the thermal character of the radiation is also dealt with. 
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O -1 : I. INTRODUCTION 

U 

Any black hole in the real universe is necessarily a dynamical one, i.e. it is neither stationary nor 
asymptotically flat. Therefore, its horizon has to be defined locally. The need to understand such 
dynamical and cosmological black holes has led to a revival of discussions on the concepts of black hole 
, itself, its singularity, horizon, and thermodynamics Indeed, the conventional definition of black holes 

implies an asymptotically flat space-time and a global definition of the event horizon. The universe, 
however, is not asymptotically flat and a global definition of the horizon is not possible. The need for 
local definition of black holes and their horizons has led to concepts such as Hayward's trapping horizon 
, and Ashtekar and Krishnan's dynamical horizon Q . It is not a trivial fact that in a specific example 
of a dynamical black hole any of these horizons may occur. In addition, we do not know if and how the 
Hawking radiation and the laws of black hole thermodynamics will apply to dynamical and cosmological 
black holes. Hence, it is of special interest to find out specific examples of dynamical black holes as a 
test bed for horizon problems. 



Within a research program, we have already found an exact solution of Einstein equations based on 
the LTB solution Q representing a dynamical black hole which is asymptotically FRW Different 
horizons and local mass definitions applied to this cosmological black hole are also reported [1, Q . In 
this paper we are looking into the question of the Hawking radiation [j| and the first law of black hole 
dynamics in such a dynamical LTB black hole . 

A first attempt to look into Hawking radiation from a 'cosmological black hole' is reported in 
There the authors consider the Einstein- Straus solution and the Sultana-Dyer one as cosmological black 
holes. The Einstein- Straus solution, however, is constructed such that it represents a 'freezed-out' 
Schwarzschild black hole within a FRW universe; it can not represent a dynamical black hole. The 
Sultana-Dyer solution is reproduced by the Schwarzschild metric through a conformal factor; it represent 
a FRW universe with a fixed black hole of a certain mass within it. Again, it can not reflect the features 
we expect from a dynamical black hole with the mass in-fall within an asymptotically FRW universe. 
Therefore, it is still an open question how the Hawking radiation and black hole dynamical laws will 
look like for a dynamical black hole within a FRW universe. We will look into this question in the case 
of the LTB cosmological black hole we have found as an exact solution of Einstein equations with mass 
in-fall, and without using a cut-and- paste technology of manifolds 

Now, Hawking's approach [8| based on quantum field theory, and applied to quasi cosmological black 
holes [9( is not a suitable method to calculate Hawking temperature in the case of proper dynamical 
black holes where one has to solve the field equations in a dynamical background. In such cases, like 
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the LTB black hole [5| , one should look for alternative approaches allowing to calculate the temperature 
of the Hawking radiation and the surface gravity. The so-called Hamilton- Jacobi (H-J) approach 11 1 
has initiated methods suitable in such cases. The method, however, suffers from not being manifestly 
covariant. Using Kodama's formalism [l9| for spherically symmetric space-times, and based on the 
Hamilton- Jacobi method, Hayward [HI has formulated a covariant form of the H-J approach and studied 
the quantum instability of a dynamical black hole [ij]; [H| • 

Another useful method is the semi-classical tunneling approach to Hawking radiation due to Parikh 
and Wilczek (PW) Ql HO, [HI 

using radial null geodesies. We will apply this radial null geodesic method, 
formulated originally for the static case, to our dynamic cosmological black hole and compare the result 
with that of the Hamilton- Jacobi method. The method, written in a gauge-invariant form, has led first 
to a Hawking temperature twice as large as the correct one (l^ |. It has, however, been shown that the 
correct Hawking temperature is regained by taking into account a contribution from the time coordinate 
upon crossing the horizon [13j . 

Among different definitions for the surface gravity of evolving horizons proposed in the past, the one 
formulated by Hayward 18| based on the Kodama's formalism [19( is the most suitable one to be used in 
our dynamic case. This leads us to the first law of black holes compatible with the Hawking temperature 
calculated by two different methods mentioned above. 

The question of thermal character of the Hawking radiation will also be discussed within this formalism. 
As Parikh and Wilczek have already pointed out, the Hawking radiation is non-thermal when the energy 



conservation is enforced [lCt \2l\. As a result it has been shown by Zhang et al 22] that the total entropy 
becomes conserved and the black hole evaporation process is unitary. 

We will introduce in section II the LTB black hole. In Section III, the covariant Hamilton- Jacobi tunneling 
method is introduced and applied to the LTB black hole. Section IV is devoted to the radial null geodesic 
method and its application to the LTB black hole. In Section V, the first law for the LTB black hole is 
derived. The thermal character of radiation is considered in section VI. We then conclude in section VII. 



II. INTRODUCING THE LTB BLACK HOLE 

The cosmological LTB black hole is defined by a cosmological spherical symmetric isotropic solution 
having an overdense mass distribution within it [f|. The overdense mass distribution collapses due to 
the dynamics of the model leading to a black hole at the center of the structure. The LTB metric is 
the simplest spherically symmetric solution of Einstein equations representing an inhomogeneous dust 
distribution It may be written in synchronous coordinates as 

ds 2 = -dt 2 + — dr 2 + R(t, r) 2 dQ 2 , (1) 

1 + f{r) 

representing a pressure-less perfect fluid satisfying 

, , 2M'(r) •, „ 2M 

'< r «'> = W' R = f+ ir- (2) 

Here dot and prime denote partial derivatives with respect to the parameters t and r respectively. The 
angular distance R, depending on the value of /, is given by 

M ... 
R = — — (1 - cosT]{r,t)), 

f-f)3/2 

r) — sin?j = — — (t — t b (r)), (3) 
i?=(-/) 1/2 ^^, (4) 

1 — COST] 

for / < 0, and 



R=(^M)i(t-t b )i, 



(5) 
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for / = 0, and 

M 

R = —(coshj](r,t) - 1), 

J3/2 

sinhrj-ri = — j-{t-t b (r)), (6) 

for / > 0. 

The metric is covariant under the rescaling r — > f (r). Therefore, one can fix one of the three free functions 
of the metric, i.e. ib(r), f(r), or M(r). The function M(r) corresponds to the Misner-Sharp mass in 
general relativity [6]. The r dependence of the bang time corresponds to a non-simultaneous big 
bang- or big-crunch-singularity. 

There are two generic singularities of this metric, where the Kretschmann scaler and Ricci one become 
infinite: the shell focusing singularity at R(t,r) = 0, and the shell crossing one at R'(t,r) — 0. However, 
there may occur that in the case of R(t, r) = the density p — -jjrw and the term remain finite. In 
this case the Kretschmann scalar remains finite and there is no shell focusing singularity. Similarly, in 
the case of vanishing R' the term may remain finite leading to a finite density and no shell crossing 
singularity either. Note that an expanding universe means generally R > 0. However, in a region around 
the center it may happen that R < 0, corresponding to the collapsing region. 

The LTB metric may also be written in a form similar to the Painleve form of the Schwarzschild metric. 
By taking the physical radius as a new coordinate using the relation dR — R'dr + Rdt one obtains 

, 9 / R? - >. , 9 dR 2 2R , „ , 
ds 2 = - I)* 3 + ^ - j+jdRdt 

+R{t,r) 2 dn 2 . (7) 

The (t, R) coordinates are usually called the physical coordinates. In the case of / = 0, the metric is 
quite similar to the Painleve metric. 



III. HAMILTON- J ACOBI METHOD 



The Hamilton- Jacobi method to calculate the Hawking radiation uses the fact that within the WKB 
approximation the tunneling probability for the classically forbidden trajectory from inside to outside 
the horizon is given by 

r oc exp f-|lm , (8) 

where S is the classical action of the (massless) particle to the leading order in h [ll|. Note that 
the exponent has to be a scalar invariant, otherwise no physical meaning could be given to T. If, in 
particular, it has the form of a thermal emission spectrum with 21m S = /3u, then both the inverse 
temperature j3 and the particle's energy ui have to be scalars; otherwise no invariant meaning could be 
given to the horizon temperature, which would then not be an observable. 

Now, let us use the Kodama formalism to introduce the necessary invariant quantities [T3|. Any spher- 
ically symmetric metric can be expressed in the form 

ds 2 =j lJ (x i )dx i dx j +R 2 {x i )dn 2 , i,je {0,1}, (9) 

where the two-dimensional metric 

d>y 2 = llj {x l )dx l dx 1 (10) 

is referred to as the normal metric, x l are associated coordinates, and R(x l ) is the areal radius considered 
as a scalar field in the normal two-dimensional space. Another relevant scalar quantity on this normal 
space is 

X {x)=j ij {x)d i Rd j R. (11) 
The dynamical trapping horizon, H , may be defined by 

X(x)\ =0, d lX \ ^0. (12) 
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The Misner-Sharp gravitational mass is then given by 

M{x) = \r(x) (1 - x (x)) , (13) 

which is an invariant quantity on the normal space. In the special case of the LTB metric this reduces 
to the Misner-Sharp mass M(r). Note also that on the horizon M\h = m = Rnj2. Now, there is 
always possible is such a spherically symmetric space time to define a preferred observer and its related 
invariant energy corresponding to the classical action of a particle. The Kodama vector (l9| . representing 
a preferred observer corresponding to the Killing vector in the static case, for the case of the LTB metric 
©is given by 

K\x) = -L= ''/),/.• K e = = K* . (14) 

Using this Kodama vector, we may introduce the invariant energy associated with a particle by means of 
its classical action being a scalar quantity on the normal space 

u = K l d t S. (15) 

Note that during the process of horizon tunneling uj is invariant independent of coordinates and is regular 
across the horizon. In the case of the eikonal approximation for massless wave field (geometric optics 
limits) , which plays an important role in calculating the Hawking radiation using tunneling method (2pj , 
the classical action S for the massless particle satisfies the Hamilton- Jacobi equation 

-•■<),Si),S 0. (16) 

The relevant imaginary part of the classical action along the 7 null curve is calculated in [li} . where it 
has been shown that the tunneling rate ([8]) is valid for the future trapped horizon, and 

Im S = Im dx l d l S^j = ^ , (17) 

where ujh is the Kodama energy and kr is the dynamical surface gravity associated with the dynamical 
horizon: 



h 2^=7 



! <•;,-( v jY J d 3 R) (18) 

H 



These are scalar quantities in the normal space. Therefore, the leading term of the tunneling rate is 
invariant, as it should be for an observable quantity. The particle production rate then takes the thermal 
form r ~ e~T with 

tlKH 



A. Application to the LTB back hole 



Assume the cosmological LTB black hole which has an infinite redsift surface satisfying the eikonal 
approximation condition for the Hamilton- Jacobi equation (|16[) . It has been shown in [3] that the apparent 
horizon of this LTB black hole in its last stages is a slowly evolving horizon with the least mass in-fall 
due to the expanding background preventing the mass infall to the central black hole. 
Now, from equations ([T]) and ([2]), and the definition of the surface gravity ((T8]) . we obtain 

mm! 1 ml 

KH = lP~ 2RW = 2R~ 2RIV 1 ( ' 

where m is the Misner-Sharp mass on the horizon. Using this expression for the surface gravity, we obtain 
the Hawking temperature according to the Hamilton- Jacobi tunneling approach: 



= ii.j.^, ',' i. (;n) 
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To relate this result to the temperature seen by the Kodama observer, we calculate first the Kodama vector 
(fT4f. It is given by K* = ^-{R\-R). The equation \K\ = = \fl + f - R 2 = ^1 - ^ 

shows that the Kodama vector is a null vector on the horizon. The corresponding velocity vector is then 
given by K l — j^. We then obtain the frequency measured by such an observer as ui — K l diS. The 
emission rate will take the thermal form 

Tcxe~f (22) 
which defines the temperature T. It is then easily seen that the temperature for this observer is 

T= ,1 , (23) 

-i 2m 

1 R 

which diverges at the horizon. The invariant redshift factor , 1 is the same factor which appear in 

\ 1 -ir 

the light frequency on the horizon showing an infinite redshift to the observer in the infinity. Then T 
itself, being finite at the horizon, may be interpreted as the redshift-renormalized temperature [l4j . 

IV. RADIAL NULL GEODESIC APPROACH 

There is another approach to the Hawking radiation as a quantum tunneling process using WKB 
approximation for radial null geodesies tunneling out from near the horizon [101 ]. The imaginary part of 
the action is defined by 



Tout rout f Pr , . , 

ImS = Im / p r dr = Im / / dp r dr 

Jr- Jr- JO 



in m 
H 



/out r" -dH' , 

Im / I — : — dr, (24) 



in 



using the Hamilton equation r — j^\ r with H being the Hamiltonian of the particle, i.e. the generator 

of the cosmic time t. Now, taking the tunneling probability as T ~ e _ 's^ m ' s , being proportional to the 
Boltzmann factor e~T , we find the Hawking temperature as 

r -£s- < 25 > 

It is easy to show that for a Schwarzschild black hole one obtain the correct expression Th = -g^jj ■ 

It has been pointed out in that 21m J^ out p r dr is not canonically invariant and thus it does 

not represent a proper observable. The object which is canonically invariant is Im § p r dr, where the 
closed path goes across the horizon and back. Using this invariant definition and T ~ e~ Im $ P rdr , the 
Hawking temperature is found to be twice the original temperature. This discrepancy in the temperature 
has been resolved by considering a temporal contribution to the tunneling amplitude. In the case of 
the Schwarzschild black hole, the temporal contribution to the action is found by changing Schwarzschild 
coordinates into Kruskal-Szekeres coordinates and then matching different Schwarzschild time coordinates 
across the horizon [13j. 

A. Application to the LTB black hole 

We use the LTB space-time in physical coordinates. In this case the outgoing and ingoing null geodesies 
are given by 

dR 



— = (R±^TTf), (26) 

where the plus sign refers to the outgoing null geodesies. Now, expanding the above equation around the 
horizon, we obtain 



« - vn7-vr+7,/i R - 2m 



dt v J v y R{l + f) 

R - 2m (27) 



G 



where R — 2m is assumed to be a non-zero small quantity. Using the results for the radial null geodesies, 
we obtain for the imaginary part of the action corresponding to the LTB metric 

f R ont 

ImS = Im / PrAR 

J R- 

m 

Im / / — ^dR 



R in Jo R 



fl in JO 

where we have used the above expansion for 4^ up to the first order of R — Rh = R — 2m. The 
corresponding Kodama invariant energy of the particle to — K l Si = \/l + fdtS is then calculated to be 
w = (VI + 7)dH'. We then use the expansion of R — 2m(r) in the form 

R-2m(r) = 

(1 _ 2 g )(fl _ H|fi) _ 2 ^r> (i _ iW . 



di 



Changing the variables (t, R) to (t, r), makes it easier to use the expression for the surface gravity in 
synchronous coordinates. Putting = ^ili = m >&. \ R=conU ^\ R=cont = -1 and = 2R in 



and (f28j) . we obtain for the imaginary part of the action 



ImS= ™ , =™, (30) 

\ 2R 2RR> I H 

leading to the Hawking temperature ()20[) : 

T = ^r- (31) 

This Hawking radiation temperature is the same as the one calculated in the previous section using 
Hamilton- Jacobi approach [15| . The definition of the surface gravity used here has been essential to 
arrive at this result indicating it to be more useful than the other definitions in the literature [l6j . 
Note that according to @ and ([7]), the term R — 2m and accordingly the metric factor of dt 2 , i.e. j^j — 1, 

vanishes on the horizon, leading to ^ = 2R. A fact which has not to be assumed while calculating the 
radial null geodesies, as has been done in Otherwise it would result in the vanishing of the imaginary 
part of the action and therefore no tunneling. 

V. FIRST LAW OF THE DYNAMICAL LTB BLACK HOLES 

Using the tunneling approach for the Hawking radiation, we formulate now a first law of the LTB black 
hole. Consider the following invariant quantity in the normal space (|33l) : 

= 7 «t w , (32) 

where TJy is the normal part of energy momentum tensor. Now, using the invariant surface gravity 
term in LTB given by (|2"0"|) , it is easy to show that on the dynamical horizon of our LTB black hole we have 

— 77T) ^R H T [ h ] , (33) 
ZRh 

where we have used Too = ~P according to the section (3.1) and Einstein equations ([2]). The horizon 
area, the areal volume associated with the horizon, and their respective differentials are then given by 

A H = 4irR 2 H , dA H = 8irR H dR H , (34) 
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Substitution from above leads to 



Vh = t^R 3 h , dV H = AirRjjdRH ■ (35) 



dV H . (36) 



Introducing the Misner-Sharp energy at the horizon, i.e. m = i?#/2, this can be recast in the form of a 
first law: 

dm = ^d(^f) - T^dV H = ^ds H T™dV H , (37) 

Z7T \ 4 J Z7T 

where sjj = Ah /47i generalizes the Bekenstein- Hawking black hole entropy. 

To see how this black hole first law is related to the Hawking radiation, we concentrate on two conserved 
currents which can be introduced in for LTB black hole. The first one is due to the Kodama vector K a , 
and the corresponding conserved charge given by the area volume V = J K a da a — 47ri? 3 /3, where da a 
is the volume form times a future directed unit normal vector of the space- like hypersurface a a . The 
second one may be defined by the energy-momentum density j a = T£K b along the Kodama vector, and 
its corresponding conserved charge E = — J j a da a being equal to the Misner-Sharp energy. The total 
energy inside the apparent horizon can then be written as Eh = R\h/2, which is the Misner-Sharp 
energy at the radius R — R\h of our LTB black hole. The energy outside the region can be expressed 
as E > h = — J T£K b d(j a , where the integration extends from the apparent horizon to infinity. We may 
therefore express the total energy of the spacetime as 

E t = R\ H /2- I T b a K b da a . (38) 

J a 

Consider now a tunneling process. The initial state before the tunneling defined by R — R\h having 
an energy E\, and the final one after the tunneling defined by R — R\h + SR\h having the energy E' t . 
According to the energy conservation, the Kodama energy change between the final and initial states of 
the tunneling process is then calculated to be 

w = dE t = E{ - E\ = 5 -^- - P dV. (39) 

Substituting from Q30p and Q36p in the above equation, we obtain the tunneling rate V ~ e f>~^ = 
eTr f»i dA H — g As^ w j^j 1 a s _ S j. _ s . being the entropy change. Our discussion shows that the tunneling 

rate arises as a natural consequence of the unified first law of thermodynamics dEn = Tds — T H dVu at 
the apparent horizon. 



VI. NON-THERMAL RADIATION FROM THE LTB BLACK HOLE 



The question of how the formulas for black hole radiation are modified due to the self-gravitation of 
the radiation is dealt with by Kraus and Wilczek in [2l[ . There it is shown that the Hawking radiation is 
non-thermal when the energy conservation is enforced. The particle in the particle hole system is treated 
as a spherical shell to have a workable model with the least degrees of freedom. The radiating black hole 
of mass M will be modeled as a shell with the energy u> around the hole having the energy M(r) — u>. 
We are going to adapt this model to our LTB black hole. Note first that due to the fact that our LTB 
black hole model is asymptotically Friedman-like, the ADM mass is not available to be fixed there. We 
therefore turn to the quasi-local Misner-Sharp mass. Let us then calculate the tunneling amplitude using 
this modification: 



ImS 



Im 



'out 



PR.dR 



= Im 



out 



H 



-dW 
R 



dR 



,2(M(r)-w) r H jtt, 

ha dR 



(40) 



Carrying out the first integral, we obtain 

" H TldH' 



Im5 f = 



( 1 m_ 'v 

UM-4F (4M-4H')ff > 

(41) 



( 4M 4Mff ) Kff 

where we have used the fact ^7 << 1 in the last step to expand the integrand and carry out the integral. 
The result shows the non-thermal character of the radiation. 

It has been shown in 0] that the boundary of the LTB black hole becomes a slowly evolving horizon for 
R' >> 1, with the surface gravity being equal to kh = jjj, and an infinite redshift for the light coming 
out of this horizon. Therefore, using the above equation, the tunneling probability has the same form as 
in the case of Schwarzschild, i.e. T ~ exp[— 8ttuj(M — £)]. Having this form of the tunneling probability 
we may refer to Zhang et al. [HJ who have shown that this form of the non-thermal radiation leads to 
the conservation of the total entropy. 



VII. CONCLUSIONS 



Within a research program to understand more in detail the LTB cosmological black hole, its 
similarities and differences to the Schwarztschild black hole, we have calculated the Hawking radiation 
from this dynamical black hole by using the tunneling methods suitable for dynamical cases. It turns out 
that for the LTB black hole the Hamilton- Jacobi and the radial null geodesic-method both lead to the 
same tunneling rate. It turns out that it is not the event horizon but the future outer tapping horizon 
that contributes to this Hawking radiation. 

Formulation of a first law for the LTB black hole and the tunneling amplitude show that the tunneling 
rate has a direct relation to the change of the LTB black hole entropy. Assuming the energy conservation 
for the LTB black hole's slowly evolving horizon, we show that the radiation is non-thermal and that 
the entropy is conserved during the radiation. 
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